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PREFACE 

Irregularities  in  the  earth’s  gravitational  field  affect  the 
motions  of  bodies  orbiting  the  earth,  The  purpose  of  this  study  is 
to  determine  the  effect  upon  the  earth'  6  gravitational  field  of  assumed 
variations  in  mass  distribution  between  oceanic  and  continental  areas. 
The  results  should  be  of  interest  to  agencies  and  persons  concerned 
with  terrestrial  structure  and  gravity,  satellite  orbits,  and  geodetic 
methods. 


V 


SUMMARY 


The  configuration  of  the  earth' s  potential  field  has  teen  de¬ 
termined  frcm  the  orbits  of  near-earth  satellites.  This  potential 
can  be  represented  by  a  series  expansion  in  spherical  harmonic s. 

The  values  of  certain  coefficients  (Jq)  for  terms  in  this  series  have 
been  calculated.  The  coefficient  reflects  the  oblateness  of  the 
earthj  the  cause  of  the  J0  coefficient  has  not  been  definitely  estab- 
li shed. 


The  potential  field  can  be  expressed  by  surface  gravity  anomalies. 
Gravity  anomalies  result  from  inhomogeneities  in  mass  distribution  with¬ 
in  a  body.  Such  differences  in  mass  distribution  exist  between  oceanic 
and  continental  areas  within  the  crust  and  tipper  mantle  over  large  areas 
of  the  earth' s  surface.  This  Memorandum  presents  the  results  of  an  in¬ 
vestigation  of  the  gravitational  fields  of  earth  models  performed  to  de¬ 
termine  the  nature  of  gravity  anomalies  to  be  expected  from  oceanic  and 
continental  areas.  Using  conventional  crustal  structure,  the  earth-model 
results  indicate  that  a  surface  anomaly  exists  which  is  positive?  over  con¬ 
tinental  bloc!;:;,  and  that  its  Magnitude  increases  with  increase  j:i  tne 
assumed  depth  to  isostatic  compensation  ana  with  increase  in  positive 
topographic  relief. 

Application  of  the  earth-mo  u:l  -anomaly  to  the  :bunl  distribution  of 
continents  and  oceans  results  in  a  ;al.Milnteu  v-J.uc  for  the  gravitational 
anomaly  which  is  of  the  sonic  order  of  rr^nituie  as,  bat  of  o; -os ;  be  sign 
to,  that  derived  fru:.’  satellite  orbits.  Conditions  that  right  re  -on-  i  1-- 
till.’-  .is  T .pan  :y  ’.re  presented. 
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LIST  OF  SYMBOLS 

a  =  distance  frcm  the  center  to  the  inner  surface  of  a  spherical 
shell 

b  =  distance  frcm  the  center  to  the  outer  surface  of  a  spherical 
shell 

G  =  gravitational  constant  (6.673  x  10~^cm\  ^sec  2) 

g  =  gravitational  attraction 

=  gravity  ancmaly  at  earth’ s  surface 

Ag  =  gravity  ancmaly  of  earth  models 

J  =  coefficient 
n 

M  =  mass  of  the  earth  (5*975  x  I027g) 

Pn  =  Legendre  polynomial 

R  =  equatorial  radius  of  the  earth  (6378.1  ±  0.1  km) 

r  =  distance  from  center  of  body  to  the  point  where  the  value  of 
potential  or  gravity  is  calculated 

t  -  thickness,  Ion 

U  ----  gravitational  potential 

“  =  colatitude 


p  -  average  density 


I.  INTRODUCTION 


The  configuration  of  the  earth1  s  potential  field,  which  relates 
directly  to  the  gravitational  field,  has  been  determined  frcsn  the 
+  psfl-t-.pnit.es.  In  this  study,  the  gravitational 

fields  of  certain  earth  models  were  investigated  to  determine  how 
these  fields  would  be  affected  by  large  areal  variations  in  mass 
distribution  within  the  earth’s  crust  and  upper  mantle.  The  results 
obtained  from  the  models  were  then  compared  with  the  satellite- derived 
data. 

The  gravitational  potential  of  the  earth  may  be  expressed  as.  an 
expansion  in  terms  of  spherical  harmonics^ ^  as 


U 


GM 

R 


Cr> 


nf  1 


P 

n 


(cos  0) 


l 


n  v  2 

where  G  is  the  gravitational  constant;  M,  the  mass  of  the  earth;  R, 
the  earth’s  equatorial  radius  (R  =  6378. 1  L  0.1  ton);  r,  the  distance 
frem  the  center  of  the  earth  to  the  point  at  which  the  potential  is 
being  calculated;  p  ,  the  Legendre  polynomial ;  and  the  colatitude. 
The  J  coefficients  determine  the  deviations  of  the  gravitational  po¬ 
tential  from  a  spherical  surface.  The  values  of  th--  J  coefficients 

n 

have  been  calculated  from  the  orbits  of  near-earth  satellites.  Oincc 
the  equatorial  plane  is  chosen  to  pass,  through  the  earth1  s  center  of 
mass,  is  zero.  Depending  on  the  satellite  data  used  and  the  method 
of  calculation,  the  value  of  J0  is  about  (1082.79  +  0.15)  >:  10"' 

c_ 

„  #  f  1 

and  J.  is  about  (-2.h  +  0.3)  x  10  .  A,c  Tin-  coefficient  ru¬ 


th- •  earth,  which  gives  a  flattening  of  about 


fleets  the  oblntrnec:: 


a 


1/298.2.  The  coefficient  reflects  a  latitudinal  asymmetry  whose 
cause  has  not  "been  definitely  established.  The  form  of  the  equation 
indicates  that  U  is  zero  at  infinity  and  that  U  increases  in  a  positive 
direction  as  r  decreases;  this  equation  is  valid  for  U  exterior  to  the 
body  whose  potential  is  being  calculated. 

If  the  earth  were  a  perfectly  spherical  body  composed  entirely 
of  concentric  spherical  shells  of  uniform  density,  the  gravitational 
attraction  of  the  earth  would  be  the  same  for  any  point  on  a  spherical 
surface  at  a  given  distance  exterior  to  the  earth.  Thus,  for  any  refer¬ 
ence  sphere  the  value  of  gravity  would  be  the  same  for  any  point  on  the 
sphere,  and  no  anomalies  would  exist.  But  on  and  below  the  surface  of 
the  earth,  areal  variations  in  the  density  of  materials  are  known  to 
exist.  In  this  case,  values  of  gravity  over  the  surface  of  the  refer¬ 
ence  sphere  vary;  the  corresponding  deviations  from  an  arbitrarily 


chosen  mean  value  are  called  gravity  anomalies. 

According  to  Munk  and  MacDonald, ^ ^  the  J  coefficients  can  also 

*  n 

be  expressed  as  gravity  anomalies  .at  the  earth’s  surface  by  the  formula 


P  (cos  0 ) 

g  =  -(n  -  1)  (G?-2?  )  J  -2 — , - 

en  n  1 1  n  ^GOS  3 Jl 


* 

One  ngl  =  0.001  gal;  one  gal  =  1  cm  sec 
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The  interior  of  the  earth  consists  of  three  major  structural 
units:  the  core,  the  mantle,  and  the  crust  (Fig.  1).  The  core,  which 
extends  from  the  center  of  the  earth  to  an  average  radius  of  approxi¬ 
mately  3^71  km,  has  an  inner  solid  zone  and  an  outer  liquid  zone.  The 
mantle  is  solid,  hut  due  to  high  temueratures  and  nrpswirpK  ■»  +  * <= 
capable  of  plastic  flow  over  long  time  periods.  It  extends  from,  the 
outer  boundary  of  the  core  to  a  radius  of  from  63OO  to  6360  km.  The 
crust,  extending  from  the  top  of  the  mantle  to  the  surface  (mean  radius, 
6371  km),  is  a  relatively  thin  and  rigid  unit. 

The  density  of  the  earth  decreases  with  distance  from  its  center. 
Major  discontinuities  occur  between  the  core,  mantle,  and  crust;  les¬ 
ser  discontinuities  are  present  at  several  other  depths.  Within  the 
core  and  lower  mantle,  the*  density -layering  is  considered  to  approximate 
concentric  spherical  surfaces.  Within  the  crust,  and  most  probably  with¬ 
in  the  upper  mantle,  the  density  pattern  is  irregular  due  to  geological 
processes.  A  well-established  and  very  distinct  variation  occurs  in  Hi'- 
density-layering  of  oceanic  and  continental  crusts.  The  oceanic  crust 
is  considerably  thinner  than  the  continental  crust,  and  it  has  a  dif¬ 
ferent  composition  (Fig.  2).  From  the  surface  downward,  the  oceanic 
crust  is  composed  of  a  rater  layer  (p  r-  1.03,  t  ~  5  km),  a  s  dJncnt 
layer  (p  ”  2.32,  t  ~  1  km),  and  a  basaltic  layer  (p  '*  3-00,  t  —  'j  km). 

The  total  thickness  of  this  crust  averages  about  11  km.  The  contia'iilal 
crust  is  composed  chiefly  of  an  upper  (granitic  layer  (p  -  2.*  :)  and  a 
lover  basaltic  layer  (p  3-00).  The  total  thickness  of  this  cru: t 
varies  considerably,  as  it  nay  range  frem  about  2C  km  under  lowland:, 
to  about  r>0  -  iC  km  under  large  mountain  r-'ing*  s.  The  averagi-  thlckn*  :: 


for  continents  is  about  33  ^ir.. 


Oceanic  crust 
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The  preservation  of  these  crustal  differences  is  explained  by 
the  theory  of  isostasy.  The  earth  is  assumed  to  be  in  isostatic  equi¬ 
librium,  vhich  means  that  at  a  "depth  of  compensation"  D  hydrostatic 
equilibrium  prevails;  thus  at  D  any  tock  unit  is  under  the  same  pressure, 
regardless  of  whether  the  unit  is  under  mountain,  lowland,  or  ocean. 

The  depth  of  compensation  is  believed  to  occur  as  deep  as,  and  probably 
below,  the  base  of  the  crust. 

The  existence  of  isostasy  in  con  finned  l»y  gravity  and  geodetic 

measurements.  An  e sample  of  the  mainton'mce  of  isostasy  occurs  during 

* 

the  formation  of  large  deltas:  The  geolo(*vlc  data  show  that  the  rock 

floor  under  the  delta  subsides  gradually  under  the  increasing  load  of 
deltaic  sediments;  however,  the  absence  of  large  anomalies  in  gravity 
observations  available  in  delta  areas  appear::  to  indicate  that  the  accre¬ 
tion  of  load  on  the  subsided  area  Jr-  Mosel y  compensated  by  the  subsidence 

(4) 

rate.  Another  e  : ample  is  the  recoil  of  the  areas  unloaded  by  the  melt¬ 

ing  of  the  Pleistocene  ice  sheets.  Uplift  motions  are  still  continuing 
in  the  Ccnndinavi eji  area  and  around  the  Great  Lakes.  Geodetic  methods 
ah ov  that  the  middle  part  oT  the  Gulf  of  Bothnia  in  Gcandinavia  has  an  up¬ 
lift  rate  of  90  cm  per  100  year:  .  The  -associated  gravity  anomaly  is  de¬ 
cidedly  negative,  indicating  that  not  all  of  the  :  ubcruntal  mass  that 

(  i|  r  ] 

flowed  outward  during  the  glacial  period  has  yrt  hud  time  to  move  back.  f  J 
The  densities  and  thicknesses  of  the  various  materials  that  com¬ 
pose  the  crust  arc  determined  by  su'-h  i-an-*iMe  data  as  surface  samples, 
drilled  cores,  and  interpretations  of  neioric  wave  velocities.  The 
degree  of  variability  of  tJ.es*-  vulu»u  is  confined  within  rather  narrow 
limits.  As  the  depth  fro;.  the  nurlV’e  increases,  i  1 1 a  ■  eroj.e  sparser 
and  less  reliable,  so  that  '-orrosponiM r+\  insert  ai  nti**B  increase. 

*do me  controversy  •  v-  ^  *  1  m  '  to  *  :  ■■  vil!.<*  thin  example 

as  a  demonstration  of  o 
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II.  METHOD  CF  INVESTIGATION 

In  this  study  we  are  interested  only  in  the  gravity  contrast  con¬ 
tributed  by  differences  in  typical  oceanic  and  continental  mass  distri¬ 
butions.  Therefore,  the  effects  on  the  earth’s  gravitational  field  due 
to  tne  earth's  rotation  and  oblateness,  the  elliptic! ty  of  the  equator, 
and  the  attraction  of  other  nearby  bodies  such  as  the  sun  and  moon  were 
excluded. 

Spherical  earth  models  incorporating  either  all-oceanic  or  all- 
continental  layering  were  constructed,  using  a  radius  of  6371  km. 
Complete  isostatic  equilibrium  was  assumed.  Density,  gravity,  and  pres¬ 
sure  values  for  the  core  and  mantle  are  given  in  Table  1.  The  oceanic 
and  continental  sections,  labeled  0  and  C,  respectively,  in  this  report 
were  based  on  typical  sections  used  in  studies  on  isostasy  and  the  in¬ 
terior  of  the  earth.  The  depth  of  compensation  was  varied  to  cover  the 
range  of  likely  possibilities  for  near-surface  isostatic  equilibrium. 

An  additional  set  of  continental  sections,  C’,  having  the  high 
average  topographic  relief  of  1  km  was  included.  (One  km,  which  is 
actually  greater  than  the  average  relief  of  continental  areas,  was 
arbitrarily  chosen  as  the  outermost  limit  of  possible  large-scale  va¬ 
riations  due  to  topography.)  Figure  3  shows  the  column  fir  representa¬ 
tion  of  the  density  values  and  distribution  of  these  sections. 

The  grovitational  attractions  of  each  of  the  resultant  models 
were  calculated  ucing  the  formula 

=  V3 « c  y  p  t&Lzjih 

rc 
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Table  1 


DENSITY  (AFTER  BULLEN),  GRAVITY,  AND  PRESSURE  IN  THE  EARTH 


(5) 


Depth 

Radius,  r 

Density,  p 

Gravity,  g 

Pressure,  p 

(ton) 

(ton) 

( g/cm3) 

(cm/sec2) 

(barsa  x  10°) 

33 

6338 

3.32 

983 

0.009 

80 

6291 

3.36 

984 

0.025 

80 

6291 

3.87 

984 

0.025 

200 

6171 

3.94 

rvQ 

0.071 

400 

5971 

4.06 

981 

0.149 

GOO 

5771 

4.18 

979 

0.230 

800 

5571 

4.30 

977 

0.313 

1000 

5371 

4.4l 

975 

0.398 

1200 

5171 

4.52 

974 

0.485 

1400 

4971 

6.63 

975 

0.574 

1600 

4771 

4.74 

977 

0. 666 

1800 

4571 

4.84 

980 

0.759 

2000 

4371 

4.94 

987 

0.855 

2200 

4171 

5.03 

996 

0.954 

2400 

3971 

5.13 

1010 

1.056 

2600 

3771 

5.22 

1026 

1.161 

2700 

3671 

5.27 

1041 

1.216 

2900 

3471 

5.57 

1068 

1.330 

2900 

3471 

9.74 

1068 

1.33 

3000 

3371 

9.90 

1043 

1.41 

3200 

3171 

10. 20 

1005 

1.62 

3400 

2971 

10.47 

960 

1.82 

3600 

2771 

10. 72 

913 

2.02 

3800 

2571 

10.99 

865 

2.21 

4000 

2371 

11. 16 

816 

2.40 

4200 

2171 

11.36 

767 

2.58 

44 00 

1971 

11. 54 

717 

2.75 

4600 

1771 

11.71 

670 

2.91 

’jSOO 

1571 

U.85 

632 

3.06' 

4982 

1389 

12.00 

593 

3.19 

5121 

1290 

15.01 

56U 

3.30 

5400 

971 

16.16 

457 

3.53 

5700 

G71 

17.07 

32 6 

3.72 

6000 

371 

17.65 

184 

3. 85 

6371 

0 

17.9 

0 

3.92 

to 


aTfco  bar 
.'*'.'7  atm. 


lr-  th»'  netcorologict, 
"/fv'  pr»-r;nuro  a*,  th- 


f  p 

z  bar  of  10  dyiy  G/an'  }  vhich  i  c  '  qua] 
center  ic  \huc  .  fC?  x  sOlP  er,c  unit*:. 
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where  g  is  the  gravitational  attraction)  p,  the  average  density;  b  and 
a,  the  distances  frcan  the  center  to  the  outer  and  inner  surfaces,  re¬ 
spectively,  of  spherical  shells;  and  r,  the  distance  from  the  center 
of  the  sphere  to  the  point  at  which  the  attraction  of  gravity  is  cal¬ 
culated  (r  >  b). 
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IH.  RESULTS  AND  DISCUSSION 

With  the  value  of  the  gravitational  attraction  for  the  oceanic 
section  being  used  as  a  reference  in  defining  the  surface  gravity 
anomalies  ( g_  =  0) ,  the  resulting  anomalies,  Ag,  for  the  continental 
sections  are  shown  in  Table  2. 


Table  2 

COOTniKMAL  GRAVITY  ANOMALIES 


CD 

Ag  (mgl) 

! 

CD 

Ag  (mgl) 

C33  i 

+3-7 

C33 

+4.7 

ch0 

+3.7 

O 

O 

+5-5 

So 

+4.2 

o 

VJ1- 

o 

+31.0 

c6o 

+4. 3 

j  c6o 

+61.0 

C70 

+4.9 

O 

O 

+89-0 

The  subscripts  of  C  indicate  the  depth  in  km  to  compensation.  For  the 
C  models,  Ag  is  positive  and  increases  slightly  as  depth  to  compensa¬ 
tion  increases.  For  the  C1  models,  Ag  is  also  positive,  but,  in  contrast 
to  the  C  models,  is  greater  in  magnitude  and  increases  at  a  greater  rate 
as  depth  of  compensation  increases.  Outward  fren  the  surface  of  the 
models,  the  gravitational  attraction  decreases  at  the  sane  rate  for  each 
set  of  models  having  the  come  depth  of  ccmpen cation. 

The  earth  models  indicate  that  the  magnitudes  of  Ag  are  comparable 
to  g.,  (4.8  mgl)  for  cantin  ntal  blocks  of  low  uveruge  relief  for  depths 
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of  compensation  near  to  and  greater  tba_n  70  km,  and  for  continental 
blocks  of  high  average  relief  at  shallower  depths  of  compensation. 

The  former  condition  more  closely  appro •ximo.tes  that  of  the  real  earth. 

The  model  values  are  expected  to  give  o*nly  rough  approximations  of  /^g  be- 

~  ^  \  -r  - -r,~„  *  -  —  -  -*  i  -l.  - j-  J  -  x 


perfect  sphere;  (3)  about  85  per  cent  o*f  tine  earth's  continental  area  has 
low  relief  (this  high  percentage  is  attributable,  in  part,  to  the  fact 
that  continental- shelf  areas  are  consU.erod  part  of  the  continental  blocks); 
and  (L)  rock  layers  in  the  earth  are  no  »t  concentric  shells  of  uniform  density. 

There  seems  to  be  a  discrepancy,  h.ovever,  in  the  sign  of  the  anomaly. 
3'ised  on  the  cra*th-model  anomalies  (Tab.  le  2),  positive  gravity  anomalies 
vjoula  o:cur  over  the  midlat  itudes  of  th-.o  Sorthern  Heir.i sphere  where  the 
continental  areas  ;ire  concentrated,  ;i<r  vever,  in  the  evaluation  of  the 
potential  formula  for  n  -  3,  is  negative  in  the  mi-ilntituder,  of  the 
Northern  Hemisphere,  causing  the  e qua po~ tent i ai  surfaces  to  deviate  in¬ 
ward.  Likewise,  the  sign  of  the  gravity  -tnomly  g,  (see  p.  p)  for  the 
same  latitudes  is  also  negative.  Thus,  whereas  the  model  results  sug¬ 
gest  a  deviation  in  the  shape  of  the  earth's  potential  roughly  resembling 
an  upsl  de-down  pear-shape,  the  third  zcrnrtl  harmonic,  as  derived  from 
satellite  orbits,  gives  a  r  ight-side-u::  pc.ur- shape  (st em  at  the  north  Pole). 
The  pear-shape  configuration  refers  tc  ^r«at  of  the  potential,  and  not  to 
the  geometrical  shape  of  the  e**irtix.  Th^ ■  following  explanations  have  been 
;onsider<vi  in  trying  to  reconcile  this  -  as  •  rfpan :y :  misinterpretation  of 
the  satellite  data;  error  in  the  struct-  ire  of  the  standard  sections;  var¬ 
iations  of  density  distribution  deeper  :ir  the  mantle  either  associated 
with,  or  in  iv:p*ui .lent  of,  the  distribution  of  ;or.tirr*r.tn;  or  in  ;k  of  i :;<j- 
static  equilibrium.  i^*;-TU::u  there  h-»s  t^Ken  .lore  :v:it  between 
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determinations  of  the  coefficient  made  by  several  persons  and  from 
various  satellite  orbits,  the  configuration  suggested  by  the  third 
zonal  harmonic  is  assumed  to  be  correct. 

Traditionally,  the  densities  for  typical  oceanic  and  continental 

UX  Vjob  vuv  MUO  Ulv*  -  .  -  UUJ^/  Ull  6QUl  tv.—  .. t'il'—  - 

the  base  of  the  oceanic  sediments.  Above  this  depth,  the  density  is 
greater  for  the  continental  section;  below  it,  the  density  for  the 
continents  is  always  less,  down  to  the  depth  of  compensation.  In  a 
spherical  body,  this  type  of  structure  ensures  a  positive  anomaly 
over  the  continental  section.  (Same  of  the  columns  in  Fig.  3  show  a 
short  interval  of  equal  density  associated  with  the  crossover. )  By 
changing  the  structure  of  the  typical  sections  so  that  the  mass  is  re¬ 
duced  and  the  center  of  gravity  is  lowered  underneath  the  continents, 
it  nay  be  possible  to  achieve  a  negative  anomaly  over  the  continents. 
This  change  can  be  accomplished  by  having  two  density  crossovers,  i.e., 
downward  frem  the  surface  the  continental- section  density  would  read: 
more  dense  than  the  oceanic- section  density,  less  dense,  and  more 
dense. 

Variations  of  density  deeper  in  the  mantle  could  account  for  the 
satellite-observed  potential  configuration.  If  these  variations  are 
independent  of  the  placement  of  continents,  then  some  neehaniim  other 
than  hydrostatic  equilibrium  is  needed  to  ex-plain  them.  If  they  are 
associated  with  continents,  density  variations  would  be  postulated  to 
occur  at  greater  depths  in  the  mantle,  perhaps  to  200  km,  or  even  deeper. 
Models  based  on  the  latter  supposition  would  be  speculative. 

If  the  requirement  that  the  earth  be  in  isostntic  equilibrium 
were  not  strictly'  adhered  to,  standard  sections  could  be  retained  by 


postulating  seme  sort  of  mass-deficiency  in  the  mantle  under  the  con¬ 
tinents.  This  might  be  accounted  for  should  crustal  rigidity  allow 
material  to  be  eroded  from  continents  and  '-.posited  into  oceans  at  a 
rate  greater  than  that  at  which  the  crust  and  upper  mantle  can  readjust 
toward  complete  isostatic  equilibrium.  Treatment  of  this  case  by  earth 
models  would  be  difficult. 


11 


IV.  CONCLUSIONS 

Layered  sequences  of  large  areal  extent  having  different  mass 
distributions  in  isostatic  equilibrium  can  cause  surface  gravity 
anomalies.  Earth  models  using  typical  oceanic  and  continental  sec¬ 
tions  result  in  a  positive  ancenaly  over  continents.  This  anomaly 
would  produce  a  potential  configuration  inverse  to  that  indicated  by 
the  third  zonal  harmonic  of  the  earth*  s  potential  field  as  derived 
from  satellite  orbits. 

Assuming  the  interpretation  based  on  satellite  orbits  to  be  cor- 
rect>  the  above  discrepancy  might  be  accounted  for  by  (l)  error  In  the 
structure  of  the  standard  sections,  (fJ)  variations  of  mass  distribution 
existing  deeper  in  the  mantle,  perhaps  on  the  order  of  100  to  ?00  km, 
or  more,  or  (3)  ieoslatic  equilibrium  sufficiently  incomplete  to  permit 
a  mass  deficiency  to  exist  under  the  continents.  The  discrepancy  seems 
to  be  most  plausibly  explained  by  variable  mass  distribution  deeper 
within  the  mantle,  \rh.h  som*  contribution  iron  incomplete  isostatic 


equi  1  ihri  tan. 
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